GMT画时间序列（三角函数）

gawk " BEGIN{for(i=0;i test.dat

gawk " BEGIN{for(i=0;i test1.dat

set ofn=1.ps

gmtset FRAME_PEN 1.6p ANOT_FONT_SIZE 19

psxy -R0/12.56/-1.5/1.5 -Ba2f0.5:"x":/a1f0.2g4:"y":WSne -W2/125/125/01ta -JX25/10 test.dat -K -Y5 > %ofn%

rem 画sin(x)曲线

gawk "{print $1,$3}" test.dat | psxy -R -W/255/0/0 -JX -K -O >> %ofn%

rem 画cos(x)曲线

gawk "{print $1,$4}" test.dat | psxy -R -W/0/255/0 -JX -K -O >> %ofn%

rem 画sin(x)+cos(x)曲线

psxy -R -JX -K -Sa0.2c -W3/255/0/0 -G255/255/0 -O test1.dat >> %ofn%

rem 以0.4为步长画cos(x)曲线上的数据点（点类型为五角星），线粗3，颜色为红色，内部充填黄色

echo 4.5 -1.2 27 0 7 LM --- y=sin(x)+cos(x)| pstext -G0/255/0 -R -JX -K -O >> %ofn%

echo 9 0.4 27 30 7 LM --- y=sin(x) | pstext -G125/125/0 -R -JX -K -O >> %ofn%

echo 1.5 0.3 27 30 7 LM --- y=cos(x) | pstext -G255/0/0 -R -JX -O >> %ofn%

del .gmt*